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ABSTRACT 

We exploit the recent observations of extremely metal-poor (EMP) stars in the Galactic halo and 
investigate the constraints on the initial mass function (IMF) of the stellar population that left these 
low-mass survivors of [Fc/H] < —2.5 and the chemical evolution that they took part in. A high-mass 
nature of IMF with the typical mass ~ 10 M@ for the stars of EMP population and the overwhelming 
contribution of low-mass members of binaries to the EMP survivors are derived from the statistics 
of carbon-enriched EMP stars with and without the enhancement of s-process elements (Komiya 
ct al. 2007, ApJ, 658, 367)). Wc first examine the analysis to confirm their results for various 
assumptions on the mass-ratio distribution function of binary members. As compared with the 
uniform distribution they used, the increase or decrease function of the mass ratio gives a higher- or 
lower-mass IMF, and a lower-mass IMF results for the independent distribution with the both members 
in the same IMF, but the derived ranges of typical mass differ less than by a factor of two and overlap 
for the extreme cases. Furthermore, we prove that the same constraints are placed on the IMF from 
the surface density of EMP stars estimated from the surveys and the chemical evolution consistent 
with the metal yields of theoretical supernova models. We then apply the derived high-mass IMF with 
the binary contribution to show that the observed mctallicity distribution function (MDF) of EMP 
stars can be reproduced not only for the shape but also for the number of EMP stars. In particular, 
the scarcity of stars below [Fe/H] ~ —4 is naturally explained in terms of the hierarchical structure 
formation, and there is no indication of significant changes in the IMF for the EMP Population. The 
present study indicates that 3 HMP /UMP stars of [Fe/H] < —4 are the primordial stars that were born 
as the low-mass members of binaries before the host clouds were polluted by their own supernovae. 
Subject headings: stars: abundances — stars: carbon — Galaxy: halo — Galaxy: formation — stars: 
mass function 



1. INTRODUCTION 

To reveal the nature of the extremely metal-poor 
(EMP) stars in the Galactic halo is the key to the under- 
standing of the formation process of the Galaxy as well 
as of the mechanism of star formation in the primordial 
and very metal-poor gas clouds. Because of the very 
low abundances of iron and other metals, these stars are 
thought to be survivors from the early days, and hence, 
are expected to carry the precious information about the 
early Universe when they were born while they reside 
in our nearby space. For a past decade, a lot of EMP 
stars have been discovered by HK sur yev dBeers et aT 



19921) and Hamburg/ESO (HES) survey (|Christlieb et al 



20011) , which enables us to use halo EMP stars as a probe 



into the early Universe. The number of known EMP stars 
exceeds several hundreds even if we limit the metallicity 
range below [Fc/H] < —2.5 and disc lose the metallicity 
distri bution function of these stars (jBeers fc Christliebl 
[20051 

One of their observed characteristics is very low fre- 
quency of stars below the metallicity [Fc/H] ~ —4. 
Despite that more than ~ 160 stars have been reg- 
istered in the mctallicity range of —4 < [Fe/H] < 
—3 by hi gh-dispersion sp ectroscopy (e.g., see SAGA 
Database: ISuda eT al. 2008), only three stars were found 
well below this metallicity; two hyper metal-poor (HMP) 
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stars of [Fe/H] < -5 HE 0107-5240 ([Fc/H] = -5.3; 
IChristlieb etall 12001) an d HE 1327-2326 ([Fe/H] = 
-5.4; iFrebel et al. 1 12001 . and one ultra metal-poor 
(UMP) star of -5 < [Fe/H] < -4, HE 0557-4840 
([Fe/H] = -4.8: lNorris et al.ll2007h . This has attracted 
wide interest, in particular, before the discovery of HE 
0557-484 in-be tween metallicity of -5 < [Fc/H] < -4. 
iKarlssonl (|2005t ) points out that such a metallicity cut- 
off can be interpreted as a result of metal spreading 
process in the sto chastic an d inho mogeneous chemical- 
enrichment model. [Karlssonj (|2006f ) then introduce a pe- 
riod of low or delayed star formation due to the nega- 
tive feedback by the Population III stars, during which 
metals spread to explain very low iron-abundance of 
HMP with the carbon yie l d from rotating s t ellar m od- 
els by iMevnet fc Maederl (|2002j ). iPrantzosI (|2003t ) ar- 
gues an early infall phase of primordial g as to alle- 
viate the paucity of low-metallicity stars. iTumlinsonl 
(|2006| ) adopts a semi-analytic approach for the hierar- 
chical structure formation and presents the model of in- 
homogeneous Galactic chemical evolution in an attempt 
of reproducing the statistical features of EMP stars and 
the re-ionization of the Universe. He addresses the con- 
straints on the IMF of population III stars, arguing high- 
mass IMF of the m ean mass at (A/) ~ 8 - 42 M Q . 
ISalvadori et al.l (j2007l ) also take a similar approach to in- 
vestigate the chemical evolution of our Galaxy with the 
mass outflow from mini-halos. In these former works, 
the low-mass star formation under the metal-deficient 
condition is introduced in rather arbitrary ways, and the 
proper explanation is yet to be devised about the nature 
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and origin of HMP/UMP stars. 

One of the decisive ingredients in studying the struc- 
ture formation and chemical evolution of Galactic halo 
is the initial mass function (IMF) of stars in the early 
days. Most of existent studies have assumed the IMF 
of EMP stars more or less similar to that of the metal- 
rich populations except for HMP and UMP stars. From 
the observations, however, we know that the EMP stars 
have the distinctive feature that a fair proportion of them 
show the surface carbon enhancement relative to iron, 
the proporti on by far larger th an the stars of younger 
populations (jRossi et al.l I1999T ) . In addition, it is re- 
vealed that the carbon-enhanced extremely metal-poor 
(CEMP) stars are divided into two sub-groups, CEMP-s 
and CEMP-nos according to the presence and absence 
of the enhancement o f s-process elements (jRvan et al.1 
I2005L lAoki et al.ll2007f ). This also forms a striking con- 
trast with the fact that their correspondences among the 
younger populations, CH stars and Ba stars, are all ob- 
served to exhibit the enhancement of s-process elements. 
Since the EMP survivors are low-mass stars, the enrich- 
ment of these elements are expected only through the 
mass transfer and/or the wind accretion from the AGB 
primaries in the binaries. Assuming this binary scenario 
and the same mechanism of c arbon enhancement as the 
stars of younger populations, iLucatello et al.l (|2005[ ) ar- 
gue an IMF with the typical mass of M mc i ~ 0.79 Af Q 
for E MP stars from the surplus of CEMP-s stars. Previ- 
ously, [Abia~e£iil- (2001) have also asserted an IMF peak- 
ing in the intermediate-mass range of 4 — 8 Mq for popu- 
lation III stars from the consideration of Galactic chem- 
ical evolution with the CN enrichment among the EMP 
stars. Furthermore, an IMF with M m d ~ 1.7 — 2.3M 
has been is discussed for the old halo stars from the 
MACHO observation in relation to the prospect that 
the observed micro-lcnsing m ay be caused by an alleged 
population of white dwarfs ([Adams fc Laughlinl [19961 : 
IChabrier et Zl[l996h . 

In order to use the carbon-enhancement to constrain 
the IMF, we should properly take into account the evo- 
lutionary peculiarity of EMP stars. It is known that for 
the stars of [Fc/H] < —2.5, there arc two mechanisms 
of carbon enhancement, while only one mechanism for 
the stars of younger populations, Pop. I and II, and also, 
that a different mode of s-process nucleosynthesis works 
(iFuiimoto et al.l 120001; [Suda et all [2004 llwamoto etHI 
12004 iNishimura et al.N2008| ). Applying these theoreti- 
cal un derstandings to the binary scenario. iKomiva et alJ 
(|2007l referred to as Paper I in the following) find that 
the IMF for EMP stars has to be high-mass with the 
typical mass of Af mc i — 10 Mq to explain the observed 
statistic features of both CEMP-s and CEMP-nos stars. 
In particular, as a consequence, it follows that the ma- 
jority of EMP stars, including CEMP stars, were born 
as the low-mass members of binary systems with the pri- 
mary stars which have shed their envelope by mass loss 
to be white dwar fs and have exploded as supernovae. 
iTumlinsonl (|2007af ) discuss the binary scenario for HMP 
stars in the similar way. 

The purpose of this paper is twofold, first to demon- 
strate the robustness of the high-mass IMF derived in 
Paper I, and then to discuss the implications to the for- 
mation and early evolution of Galaxy. In the follow- 
ing, we make a distinction between the total assembly 



of EMP stars that were born in the early Galaxy, in- 
cluding massive stars which were already exploded as 
supernovae, and the low-mass EMP stars that are still 
alive in the nuclear burning stages by calling the for- 
mer "EMP population" and the latter "EMP survivors" . 
In deriving the constraints on the IMF of stars for the 
EMP population, one has to make the assumptions on 
the binary characteristics, among which the most crucial 
is the distribution function of mass ratio between the pri- 
mary and secondary stars in binaries. Paper I adopts a 
flat distribution for simplicity. It seems plausible from 
the obse rvations of the stellar systems of younger pop - 
ulations ([Duquennov fe Mavorlll99lHMavor et al.lll992h . 
and yet, it is true that the mass-ratio distribution is yet 
to be properly established both observationally and the- 
oretically even for the binaries of younger populations. 
Several different mechanisms have been proposed for the 
binary formation, such as the fragmentation during the 
collapse and the capture of formed stars, and are thought 
to give different mas s-ratio distributions (see also e.g., 
iGoodwin et al.l [20071 and the references therein). The 
distribution may increase or decrease with the mass- 
ratio, or the two stars may form in the same IMF as 
suggested for the capture origin. In this paper, we exam- 
ine the dependence of the resultant IMF on the assumed 
mass-ratio distributions of various functional forms, in- 
cluding the independent coupling of the both stars in the 
same IMF to demonstrate that the high-mass nature of 
IMF of EMP population is essentially unaltered. 

The recent large-scaled surveys of EMP stars provide 
the additional information on the early history of Galac- 
tic halo. A fairly large number of known metal-poor 
stars (144 and 234 stars of [Fc/H] < -3 by the HK 
and HES surveys, respectively) makes it fea sible to dis- 
cuss the metallicity distribution function ([Beers et al.l 
|2005[) . Moreover, the significant coverage of celestial 
sphere (6900 and 8225 deg 2 by the HK and HES surveys 
respectively; IChristliebl 120031 : IBeers fe Christhebl l2005h 
allows to consider the total number of EMP survivors in 
the Galactic halo. We demonstrate that the latter also 
places an independent constraint on the IMF of EMP 
population in combination with the metal yields pro- 
duced by the EMP supernovae if the binary contribution 
is properly taken into account. 

We then apply the IMF, thud derived, to discuss the 
chemical evolution in which the stars of EMP popula- 
tion take part. It is shown that the resultant IMFs can 
reproduce the number and slope of observed metallic- 
ity distribution functions (MDF) for EMP stars, and 
also, to give an explanation to the scarcity and origin of 
HMP/UMP stars with the effects of hierarchical struc- 
ture formation process included. In this paper, we deal 
only with the iron production by SNe since we are inter- 
ested in the MDF, and discuss the basic characteristics 
of hierarchical structure formation by usin g simple an- 
alytic approximations. ITumlinsonl (|2007b[ ) studies the 
low-mass star formation taking into account the contri- 
bution of binary stars, but his approach is different from 
ours in the uses the hypothesized IMF with the effect of 
cosmic microwave background. In addition, he consid- 
ered only the CEMP-s stars, but not CEMP-nos stars or 
MDF. 

This paper is organized as follows. In §2, we discuss 
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the constraints on the IMF of EMP population from the 
statistics of CEMP stars and from the total number of 
EMP survivors in our Galaxies. In §3, we investigate the 
metallicity distribution of EMP stars in Galactic halo 
with the formation process of the Galaxy taken into ac- 
count. Then our conclusions follow with discussion of the 
origin of observed MDF and also of HMP stars. In Ap- 
pendix, we re-discuss the relationship between the num- 
ber of EMP survivors, estimated from the surveys, and 
the metal production by the EMP supernovae with the 
binary contribution taken into account, to demonstrate 
that they entail the same IMFs as drawn independently 
from the statistics of CEMP stars. 

2. CONSTRAINTS ON IMF OF EMP STARS 

In this section, we revisit the problem of constrain- 
ing the IMF for the stars of EMP population from the 
observations of EMP survivors, studied in Paper I. The 
method is based on the analysis of statistics of CEMP 
stars in the framework of binary scenario, and hence, in- 
volves the assumptions of EMP binary systems. We start 
with reviewing the method and assumptions used in Pa- 
per I in deriving the constraints on the IMF of EMP 
population stars. We first investigate the dependence 
of resultant IMF on these assumptions, in particular of 
the mass-ratio distribution of binary members. We then 
discuss the iron production by EMP population stars in 
relation to the total number of EMP survivors, estimated 
from the HK and HES surveys, to assess the constraints 
on the IMF through the chemical evolution of Galactic 
halo. 

2.1. Method and basic assumptions 

We give the outline of our method in studying the 
statistics of CEMP stars and chemical evolution of Galac- 
tic halo with the discussion of the assumptions involved, 
and a brief summary of the observational facts that our 
study rely on. 

2.1.1. Statistics of CEMP stars 

Our method is founded on the results of stellar evo- 
lution that the stars of [Fc/H] < —2.5 and of mass 
< 3.5 Mq undergo hydrogen mixing into the helium con- 
vection during the helium core or shell-flashes, differently 
from the stars of younger populations, Pop. I and II 
(jFuiimoto et al.l I1990L l2000l ). This triggers the helium- 
flash driven deep mixing (He-FDDM) to carry out car- 
bon to the sur face (|Hollowell et all fl990h . It is neces- 
sarily accompanied with the s-process nucleosynthesis 
in the helium convect ion as mixed protons are trans- 
formed into neutrons JSuda et ail |2004| ; llwamoto et al.l 
12004 iNishimura et al.ll2008f l. For EMP stars, He-FDDM 
works as the mechanism to enrich both carbon and s- 
process elements in their surface in addition to the third 
dredge-up (TDU), the latter of which works in the stars 
of M >~ 1.5 Mq in common with the stars of younger 
populations. 

Consequently, the origins of two sub-groups of CEMP 
stars are identified with these two mechanisms. The 
CEMP-s and CEMP-nos stars stem from the low-mass 
members of EMP binaries with the primaries in the mass 
ranges of 0.8 Mq < M < 3.5 Mq and 3.5 Mq < M < 
M up , respectively. Here M up is the upper limit to ini- 
tial mass of stars for the formation of white dwarfs. We 



take M up = 6.0M a (|Cassisi &: Castellanilll993l see also 
Siess 2007), which is also taken to be the lower mass 
limit to the stars that explode as supernova. This is the 
fundamental premise of our study. Among the CEMP- 
nos stars, there are the stars that show different charac- 
teristics such as CS22892-052 with a large enhancement 
of r-process elements. They may have different origins 
according to the scenarios suc h as proposed in connec- 
tion t o supernovae yields (e.g ., [Tsujimoto & Shi gevamal 
[20011 lUmeda fc Nomotol (20021 IWanaio et al.ll200fif ). Ac- 
cordingly, the CEMP-nos stars may be the admixture 
of the stars of different origins. Since the ratio between 
the CEMP-s and CEMP-nos stars depends strongly on 
the IMF, as shown from Paper I, however, our results 
will not affected as long as the CEMP-nos stars contains 
those with the AGB mass transfer origin that we pro- 
pose. 

For the formation of CEMP stars in the binary sys- 
tems, the initial separation, A, has to be large enough 
to allow the primary stars to evolve through the AGB 
stage without suffering from the Roche lobe overflow, 
but small enough for the secondary stars to accrete 
a sufficient mass of the wind to pollute their surface 
with the envelope matter processed and ejected by the 
AGB companion. The lower bound, A m \- a (mi,m2) 1 
to the initial separation is estimated from the stellar 
radii o f EMP stars taken from the evolutionary calcu- 
lation (|Suda fc Fuiimotol l2007f) . where mi and m-i are 
the masses of primary and secondary stars. The AGB 
star is assumed to eject the carbon enhanced matter of 
[C/H] = with the wind velocity w W i n d = 20km s _1 until 
it becomes a white dwarf, and we define CEMP stars as 
[C/Fe] < 0.5. The upper bound, A max (mi, 7712), is esti- 
mated by the amount of accreted matter calculated by 
applying the Bondi-Hoyle accretion rate, 

dt ~ A(t) 2 v Ie i(t)* v wind dt ' [ ' 

in the spherically symmetric wind from the companion, 
and v Te \ is the relative velocity of the secondary star to 
the wind. Accreted matter is mixed in surface convec- 
tion of depth 0.35 Mq and 0.0035 Mq in mass for giants 
and dwarfs, respectively. E.g., for the stellar metallic- 
ity [Fe/H] = —3.5, the mass of accreted matter has to 
be larger than 3.5 x 10" 4 M Q and 3.5 x 10~ 6 M Q , and 
hence, the upper bounds are ~ 100AU ~ 1000AU for 
dwarfs, respe ctively. It is poi nted out that the molecu- 
lar diffu sion dWeiss et al.ll200"ol ) and/or the thermohalinc 
mixing (|Stancliffe et al.l l2007| ) work in the envelope of 
EMP dwarfs to lower the surface abundance of accreted 
matter by nearly an order of magnitude. If these ef- 
fects are included, it demands a larger accreted mass, 
and hence, a decrease of upper bound of binary sepa- 
ration to the carbon enrichment for the dwarf low-mass 
members with carbon by ~ 0.5dex. This will not affect 
our results so much since the upper bound is itself suf- 
ficiently large to exceed the separation at the peak of 
period distribution (see below). In addition, it has little 
effects on the giants which we mainly deal with in the 
following because of by far deeper surface convection. 

If we specify the initial mass function, £(,(to), and 
the distributions of binary parameters, therefore, we can 
evaluate the frequency of CEMP-s and CEMP-nos stars, 
and through the comparison with the observations, we 
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may impose the constraints on the IMF and on the bi- 
nary parameters. The numbers of CEMP-s and CEMP- 
nos stars currently observable in flux-limited samples are 
given by 

M.SM @ 

V'cemp-s — hi dm 2 N s (L[m 2 ]) 

J 0.08 M G 

x / rfmi^(mi) 

J0.8M & 
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< I dmi£ b (mi) 

'3.5M 



0.08 M s 
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where h is the binary fraction: n(q) is the distribution of 
the mass-ratio, q = m 2 jm\, and f(P) is the distribution 
of the period of binaries: and N(L) is the probability of 
the stars in the Galactic halo with the luminosity L in 
the survey volume of HES survey. Note that all of them 
stem from the low- mass members of binary since it is only 
a very small fraction of stars that have experienced He- 
FDDM to develop CEMP-s characteristics in themselves 
and now stay on AGB. Similarly the total number of 
EMP survivors is given by 

fO.8M 

f/Wv = / dmN s (L[m])[(l - h)£s(m) + Mb(m)] 

J0.08M Q 

r0.8M e ,oo / n 

+f b dm 2 N s (L[m 2 }) dmi&(mi)-^, (4) 

J0.08M e J0.8M Q m i 

with the contribution of the stars born as single under 
the initial mass function £ s ). The rest of the terms give 
the number of EMP survivors formed as binary, ^binary 
The stellar luminosity and lifetime a re taken from the 
evolut ion calculation of EMP stars by iSuda fe Fuiimotol 
(|2007f ). The AGB primaries dredge up to increase their 
surface helium abundances, and hence, may cause the 
surface enrichment of helium to the companion stars in 
the binaries at the same time with the carbon enrich- 
ment, though both suffering the dilution in the envelope 
convection. The surface enrichment increases the lumi- 
nosity during the RGB evolution to shorten their RGB 
lifetime of polluted EMP stars nearly in inverse propor- 
tion of the luminosity, but the survey volume increases 
with the luminosity. For flux limited sample, the ob- 
served number of EMP giants may rather increases with 
the surface enrichment in proportion to a half power un- 
der a constant density distribution. On the other hand, 
the increase in the luminosity occurs only after the hy- 
drogen burning shell comes to take place in the shell to 
which the pollutants are carried in by the surface con- 
vection, and hence, in the later stages for metal-poorer 
stars. Accordingly, the effect of helium enhancement 
will little affect our results of constraints on the IMF 
since HES survey is thought to reach far enough that the 
spatial distribution of halo stars decreases. 

2.1.2. Model Parameters 

In this paper, we assume that the binary primary 
stars and single stars are born under the same IMF, i.e., 
£(m) = £h(m) = £,(m). For the form of IMF, we may 



well assume a lognormal function with the medium mass, 
M m( j, and the dispersion, Am, as parameters 



£ (m) oc — exp 
m 



(logm - logM md ) 5 



9A 2 



(5) 



In addition, we assume the binary fraction h = 0-5 in 
this paper. Our results are little affected by the as- 
sumption about h since not only the CEMP stars but 
also most of the EMP survivors come from the sec- 
ondary companions of binaries unless Af m( j < 0.8 Mq, 
as seen later. As for the binary period, we may 
adopt the distribution deri ved for the nearby stars by 
iDuquennov fc Mavorl (|199lh . 



f{P) oc p exp 



(logP 



2 x 2.3 2 
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where P is the period in units of days. The binary frac- 
tions and period distributions of halo stars are observed 
to be not significantly diffe r ent from those of ne arby disk 
stars ([Latham et all l2002t ICarnev et all l2003f ). Addi- 
tionally, it is shown in Paper I that this period distribu- 
tion is consistent with the observations of CEMP stars 
for the periods of P < 10 yr confirmed to date (see Fig. 3 
in Paper I). 

The mass ratio distribution is an essential factor in 
discussing the evolution of binary systems, and yet, it 
is not well understood. The mass ratio distribution 
of metal- poor halo stars is investigated observationally 
(e.g., see lGoldberg et al.ll2003l lAbtll2008fl . and yet, sub- 
ject to large uncertainties. Especially for the binary with 
intermediate-mass or massive primary stars, it is hard to 
know the mass ratio distribution from the observations. 
Theoretically, neither the fragmentation of gas cloud nor 
the accretion process onto proto-binaries are yet well un- 
derstood even for Population I stars (e .g. iBate fc Bonnelll 
Il997t lOchi et all [20051 : lMachida|[200l . In order to test 
the assumption on the mass-ratio distribution, we inves- 
tigate the constraints on the IMF for different mass-ratio 
distributions. In Paper I, the simplest flat distribution is 
assumed in Paper I among the possible distributions. In 
this paper, we test some other assumptions and discuss 
the dependence of IMF parameters on the mass ratio 
distributions, as stated in £12.21 

We may define the coupling mass distribution function, 
x(mi, m 2 ), as the fraction of the binaries with a primary 
and secondary star in the mass range of [mi, mi + dmi] 
and [m,2, m 2 + dm 2 ] (mi > m 2 ) and write it in the form; 

%(mi, m 2 )dmidm 2 = £(mi)n(q)dqdmi 

= ((mi)n(ffl2, mi)/midmidm 2 .(7) 

Here the initial mass function, £, of the primary star is 
assumed to be the same as IMF of single stars: And n(q) 
is the mass ratio distribution, for which we assume both 
extremities of increase and decrease functional forms in 
addition to the constant one, adopted in Paper I; 



1/(1 - 0.08 M Q /mi) (Case A) 
n{q) = { 2q/[l - (0.08 M Q /mi) 2 ] (Case B) 
g-7m(mi/0.08M Q ) (Case C). 



(8) 



Furthermore, we take up a different type of mass-ratio 
distribution that the primary and secondary stars in- 
dependently obey the same IMF such as assumed by 
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iLucatello et al.l (|2005f ). In this case, the coupling mass 
distribution function is given as a product of the same 
IMF as; 

x(mi, m2)dm\dm2 = 2£(mi)£(m2)dmidm,2 (Case D). 

(9) 

We shall refer to this distribution function as "inde- 
pendent" coupling. From the comparison with eq. (O, 
we may write the mass-ratio function in the form 
n(i7i2,mi) = 2m 1 £(m 2 ); it is should be noted, however, 
that the frequency of binaries with a primary star of mass 
mi is not normalized and increases with mi from zero to 
2, as given by the integral J om Me/mi n(m 1 ,m 2 )dq = 

2 J™08M r J( m ^ dm 2- 

With these specification and with the assumed mass- 
ratio distribution function, we may compute the fractions 
of EMP survivors, V'surv(M m d, Am) and of both cemps 
stars, ipCEMP-s(M md , A M ) and ^cEMP-nos(Mnd, A M ), 
and search the ranges of the IMF parameters, medium 
mass Mmd and dispersion Am, that can reproduce the 
statistics of CEMP stars consistent with observations 

2.1.3. Total iron yield of EMP supernovae 

We can pose another constraint from the total iron 
yield, Mp c , of EMP population and the total number, 
-^Vemp.g, of the giant EMP survivors. For A^emp,g, es- 
timated from the results of existent surveys, the total 
stellar mass, Memp, of EMP population for an assumed 
IMF is given by, 

M EM p(M md , A M ) = m A EM p,g//g, (10) 

where /g is the fraction of giant EMP survivors in all 
the stellar systems, born as EMP population, and m is 
the averaged mass of EMP population stars: 



With these evaluations and the observed number of 
EMP giants, we can give the total iron yield of stars 
of EMP population as a function of IMF parameters. 
The comparison with the total amount of iron estimated 
from the chemical evolution of Galactic halo may impose 
constraint on the IMF parameters. 

2.1.4. Observational constraints 

The first constraint is the number fraction of CEMP- 
s stars. The HK and HES observations tell that the 
CEMP stars wi th [C/Fe] > 1 account for 20 ~ 25% of 
EMP s tars fe.g..lBeerslll999tlR"ossi et al.|[l999HChristliebl 
I2003D . iCohen et al.l (|2005f) suggest a slightly lower frac- 
tion of 14.4% ± 4% with the er rors in the abundance 
analysis taken into account while ILucatello et al.l (|2006f ) 
obtain a larger frequency of 21% ± 2% for the HERES 
(HES r-process enhanced star) survey sample, both for 
[Fe/H] < -2.0. It is claimed that the frequency of CEMP 
is higher at lower metallicity of [Fe/H] < —2.5 but we 
have to subtract contribution from the CEMP-nos stars. 
In this paper, we adopt the observational constraint on 
the fraction of the CEMP-s stars at 10 - 25%; 



0.1 < 



V'CEMP-s(Afmd, Am) 



< 0.25. 
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The first terms of both equations denote the contribu- 
tions by the stars born as the single stars and as the pri- 
mary stars in the binaries and the second terms denote 
the contributions by the stars born as the secondary stars 
in the binaries. The mass and mass range of EMP stars 
now on the giant branch are taken to be M = 0.8 M Q 
and AMq — 0.01 M Q , based on the stellar evolution cal- 
culation of stars with [Fe/H] = —3, as in paper I. 

The massive stars of EMP population have exploded 
as supernovae to enrich the interstellar gas with metals. 
The amount of iron, Mpe.EMP, ejected by all the super- 
novae of EMP population of the total mass, Memp, is 
given by 

MfcEMP = ^^/sn(Ifc) = N EMP , G ^{Y Fc ), (13) 

m Jg 
where /sn is the fraction of the stars that have exploded 
as supernovae and given by, 



/sn = 



dmi£(m 



)[1 + A f 

mi J Mu 



n(q)dm 2 ] : (14) 



and (Yfc) is the averaged iron yield per supernova, taken 
to be (Yfv) = 0.07 Mm in the following calculations. 



VWrv(M,nd, Am) 
The second constraint is the number ratio between 
CEMP-nos and CEMP-s stars. The observed frequency 
of CEMP-nos to CEMP-s stars is ~ 1/3 or more , (e.g. , 
iRvan et al.l [20051 : lAoki et alJl200l . lAoki et al.l ((200% 
point out that it increases for lower metallicity, report- 
ing the ratio as large as 9/14 for [Fe/H] < —2.5. In 
addition, EMP stars enriched with nitrogen are found 
in number compara ble with, or more than, CEMP-nos 
stars ("mixed" stars: ISpite et al.ll2005l ). whose origin can 
be interpreted in terms of the same mechanism but with 
more massive primary companions that experience the 
AM(hot bottom burning (HBB) in the envelope of the AGB. 
Some other scenari o s for CEMP stars have been p roposed 
lUmeda fc Nomotol (|2005l) : iMevnet et all (|2006f ) but we 
assume all CEMP stars are formed in binaries with AGB 
in this paper. We adopt the observational constraint on 
the relative frequency of CEMP-nos to CEMP-s stars at 
1/3-1; 

, ,„ , V'CEMP-nos(M md , A M ) , n „s 
I/ 3 < — 7TT T-^- < 1- (16) 

We note that the above two constraints are not depen- 
dent on the total mass nor on the spatial distribution 
of the stellar halo because they are concerned with the 
relative number ratios. 

The third constraint is the total iron yield from the 
EMP population. The HES survey obtained 234 stars 
of [Fe/H] < -3 (|Beers et al.l 120051) as a result of the 
medium-resolution, follow-up observations of 40% of the 
candidates, selected by the objective-prism survey of the 
nominal area S = 8225deg 2 Jb cers fc Christliebll2005h . 
Taking the relative frequency between the giants and 
dwarfs (1 : 0.93) and the ratio of the stars of [Fc/H] < —3 
and [Fc/H] < -2.5 (6% : 20%) from their Table 3, we 
may estimate the total number of EMP giants in the 
Galactic halo in the flux limited sample at: 

^emp.g — 410 sr ^. (17) 

We assume that all giant stars in the survey areas are 
observed because of the fairly large limiting magnitude 
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of HES survey (B = 17.5), about two magnitude deeper 
than for the HK survey, and neglect the spatial distri- 
bution of EMP giant for simplicity since the sufficient 
information is not yet available (see §7.3 in Paper I for 
the detail). Then we have the total number of EMP gi- 
ants -/Vemp.g = 5.2 x 10 3 in the Galaxy. 

On the other hand, the amount of iron necessary 
to promote the chemical evolution of the whole gas in 
Galaxy of mass, M h = 10 11 M Q , up to [Fe/H] = -2.5 is 
as much as 

A/Fchaio = M^f^qKT 2 - 5 = 10 5 - 5 M Q , (18) 

and the supernovae of EMP population should have pro- 
vided this amount of iron unless there were other pop- 
ulation^) of stars which made iron without producing 
low- mass stars. Using eq. (p~3|) . this is transferred into a 
constraint on the IMF as 

M Fc . EMP = 0.07Af©x5.2xl0 3 / f/ff md ^ M > ~ 10 5 ' 5 Mq 

(19) 

The estimated number of EMP survivors may be subject 
to significant uncertainties. If we take into account the 
EMP stars in the outer halo and in the Galactic bulge 
that HES survey cannot reach, ./Vemp.g can be larger, 
which demands a smaller amount of iron produced per 
EMP survivor, and hence, a smaller number of super- 
nova, leading to a lower-mass IMF. A lower-mass IMF 
also results if the binary fraction is smaller and/or if there 
is other source(s) of iron that does not accompany the 
low-mass star formation. On the other hand, if a part 
of the supernovae ejecta is dispelled and lost from the 
Galaxy, it demands a larger amount of iron, AfF e) EMP, 
and hence, a higher-mass IMF. Despite such uncertain- 
ties both of the observations and the theoretical assump- 
tions, the constraints on the IMF derived from eq. (|19|) 
are rather robust since the ratio, /sn//gj is a rapidly 
varying function of IMF. 

2.2. Dependence on Mass-Ratio Distributions 

For the four mass-ratio distributions, formulated in 
^2.1.21 we can figure out, as the function of M m( j and Am, 
the portion of stars that survive to date (M < 0.8 M©), 
and then, the fractions of stars in these EMP survivors 
that evolved to CEMP-s and CEMP-nos stars according 
to the masses of primary stars and to the orbital separa- 
tions. Figure [1] compares the fractions of CEMP-s stars 
in the EMP survivors [T/'cEMP-s/V'surv] and the ratios of 
CEMP-nos to CEMP-s stars [f/'CEMP-nos/V'CEMP-s], pre- 
dicted with use of these four different mass-ratio distri- 
butions, as a function of medium mass M mc i of IMFs 
with the dispersion of Am = 0.33, taken to be same 
as th e present-day IMF of Galactic spheroid compo- 
nent dChabrierll200l . Fi gures [5] and [3] present the con- 
tour maps on the M m( j-AM diagram for the fractions of 
CEMP-s stars and the ratio between the CEMP-nos and 
CEMP-s stars, respectively. 

Left top panels on these figures show the results for 
the flat mass-ratio distribution of Case A, which repro- 
duces the results obtained in Paper I. In Fig. Q] the 
CEMP-s fraction peaks at M md = 4.8 M©, slightly above 
the upper mass limit of the primary stars for CEMP- 
s. Note that when the secondary mass is specified, the 
mass distribution of primary stars peaks at mass smaller 
than M m H for this mass-ratio function foe f(mi)/mi, see 



Fig. 12 in Paper I). Two ranges of M md , 0.61 - 2.8 M© 
and 7.6 — 15.3 M© (light shaded parts) gives the IMFs 
compatible with the observations, separated by the over- 
production of CEMP-s stars. The relative frequency of 
CEMP-nos to CEMP-s stars is a steep increase func- 
tion of M m( j, and excludes the lower range of M m( j 
compatible with the CEMP-s fraction. The IMFs with 
Mmd = 4.8 — 11.6 M© (dark shaded part) gives com- 
patible ratio with the observations This range of M m( j 
lies in the mass range of primary stars of CEMP-nos 
stars or even larger. Accordingly, the intersection of 
the light and dark shaded parts designates the ranges 
(M m d = 7.6— 11.6 Mq) that can explain the both statis- 
tics of CEMP stars, and hence, high-mass IMFs results 
for a dispersion Am = 0.33. 

On the M m d — Am diagram of Fig(3J the parame- 
ter space compatible with the observed CEMP-s frac- 
tion separates into two ranges for the dispersion smaller 
than Am — 0.43, converging to the narrow ranges around 
M m d — 1 and 4Af©, respectively, as Am decreases. For 
larger dispersion, on the other hand, it merges into one 
part to cover wider range. As for the ratio between the 
CEMP-nos and CEMP-s stars, Fig. [3] shows that the 
medium mass compatible with observed ratio increases 
with the dispersion to cover wider range, from Af,„d = 
3.2- 3.7 M© at A M = 0.1 through M md = 12.3- 100 Af© 
at Am = 0.54. Accordingly, for the IMFs that satisfy the 
both statistical constraints, the medium mass increases 
with the dispersion from M m d > 5.5 M© for Am = 0.22 
and beyond A/ md = 100 Af© for A M > 0.62. 

For the mass-ratio distribution function increasing 
with q of Case B (right top panel), the portion of bi- 
naries that have the secondary stars surviving to date 
decreases with the mass of primary stars in propor- 
tion to (7712/7711 )~ 2 , more steeply than in proportion to 
(7712/7711) -1 for a flat mass-ratio distribution in Case A. 
Since the average mass of the primary stars is smaller 
for a given EMP star, therefore, the fraction of CEMP- 
s stars is larger for a given M m d, and the peak shifts to 
larger M m d, as compared with Case A. In Fig.[TJ the Af m d 
of IMFs compatible with the observed fraction of CEMP- 
s stars separates into two ranges, as in Case A, but the 
in-between gap is larger; the higher mass range shifts up- 
ward in mass to greater extent (M m d = 14.5 — 28A^f©) 
than the smaller mass range shifts downward (Af m d = 
0.38-2.4 Mq). This also causes a smaller ratio of CEMP- 
nos to CEMP-s stars for a given M m d, and hence, the 
IMFs compatible with the observed ratio shift to a larger 
mass of A/ m d = 8.9 — 21.6 A/©, as compared with that for 
Case A. As a result, the IMFs consistent with the both 
statistics of CEMP stars turns out to be higher mass by a 
factor of - 2 than for Case A (M md = 14.5 - 21.6 M© for 
Am = 0.33). In the Fig. ® we see that the range of M m d, 
compatible with the observed fractions of CEMP-s star 
(shaded area), separates into two and the higher range 
shifts to larger mass for a given Am- Similarly, in the 
Fig. [3j the observed ratio of the CEMP-nos to CEMP-s 
stars also demands larger Af m d, and the range of M m d of 
IMFs compatible with the observation increases rapidly 
with Am to exceed 100 Af© for Am > 0.59. In order to 
satisfy the both conditions of CEMP star observations, 
the IMFs fall in the range of higher medium mass and in 
a rather narrow range of dispersion, lying in the param- 
eter space of Af m d > 7.1 Af©, larger by a factor of ~ 1.3 
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Fig. 1. — Constraints on the IMF of EMP population for four cases of assumption for mass-ratio distribution function. Top Ze/fCase 
A: n(q) = const., Top rightCasc B. n(q) oc q, Bottom Ze/iCase C. n(q) oc 1/q, and Bottom rightCasc D. independent coupling. Thin and 
thick solid lines denote the fraction of CEMP-s in the EMP survivors as the function of the medium mass M mt j for a fixed dispersion 
of Am = 0.33 for the EMP stars born as binaries (i/'CEMP-s (M mt j, 0.33)/i/'bmary (Mmdi 0.33)) and for the total EMP stars including the 
single stars born in an equal number to the binaries(»/!cEMP-s (A^mdi 0.33)/i/'surv(A , / m( ji 0.33)), respectively. Broken line denotes the ratio 
of CEMP-nos to CEMP-s stars (^CEMP-nos(-^mdi 0.33)/i/)cEMP-s(-^mdi 0.33)). Light and dark shaded areas denote the parameter ranges 
for the IMFs that can give rise to the observed fraction of CEMP-s stars in the EMP survivors (10 — 25%) and the observed ratio of 
CEMP-nos to CEMP-s stars (1/3 - 1), respectively. 



than for Case A, and A M > 0.22 and of A M = 0.45-0.58 
for M md = 100 M Q . 

For a mass ratio function decreasing with q of Case C, 
we see the opposite tendency of Case B (bottom left 
panels of Fig. [T] - |3J). The portion of EMP bina- 
ries whose low-mass members survive to date depends 
only weakly on the primary mass (oc log mi) so that 
the fraction of CEMP-s stars reduces because of larger 
contributions from the binaries with more massive pri- 
maries. As seen in Fig. [T] the fraction of CEMP-s stars 
in the total EMP survivors is well below the upper bound 
of the observations, and hence, the M m d compatible 
with the observations merges into one narrower range 
of M md = 1.1 - 15.6 M for A M = 0.33. The observed 
ratio of CEMP-nos to CEMP-s stars can be reproduced 
also by the IMFs with a smaller M m d by a factor of ~ 2 
than in Case A (M md = 3.2 - 7.5 Mq). Accordingly, 
the M m d for the IMFs consistent with the both CEMP 
star statistics are smaller by a factor of 1.5 — 2.4 than 
for Case A (the mass range M m( j = 3.2 — 7.5 Mq for 
A M = 0.33). In the Fig. [2J the range of M md for the 
IMFs, compatible with the observed CEMP-s fraction 
varies only little with Am, and is restricted in the range 
between M md = 1.1—23 Mq, though it separates into two 
for small Am < 0.31. As shown in Fig.[3]the dependence 
of the ratio of CEMP-nos to CEMP-s stars on Am is also 
weaker than for Case A. Consequently, the IMFs can re- 



produce the both CEMP star statistics with the mass as 
small as M m d = 3.3 Mq, smaller by a factor of ~ 0.6 than 
for Case A, but differently from the above two cases, an 
upper bound is placed at M m d = 23 Mq, regardless of 
the dispersion with a lower bound of Am > 0.21. 

Bottom right panels depict the results for the "inde- 
pendent" coupling of Case D. For this case, the number 
of EMP survivors produced per binary is independent 
of the mass, mi, of primary stars, while the binary fre- 
quency itself increases with mi. The former is similarly 
to Case C, and then, the production of EMP survivors 
from the binaries with massive primary poses a severe 
constraint on the high-mass side of IMFs. On the other 
hand, the latter favors the production of CEMP-s stars 
as compared with the low-mass binaries of mi < 0.8 Mq. 
The both shift the IMFs, compatible with the observed 
fraction of CEMP-s stars, to smaller M m d. In addi- 
tion, the single stars, born in the same number of bi- 
naries, contribute to significant fraction of EMP sur- 
vivors, increasing from 1/3 up to 1/2 for smaller M m d for 
M m d < 0.8 Mq since the low-mass binaries are counted 
as one object. As a result, the maximum fraction of 
CEMP-s stars remains below the upper limit of the ob- 
served range, which makes the M md for the IMFs that 
can reproduce the observation lie in a single range within 
a relatively small upper bound. The observed ratio of 
CEMP-nos to CEMP-s stars demands also lower-mass 
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Fig. 2. — Dependence of the fraction of CEMP-s stars on the medium mass, M m< j, and dispersion, Am of IMF of EMP population. Solid 
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numerals designate the fractions. Shaded area denote the parameter ranges for the IMFs that can give rise to the observed fraction of 
CEMP-s stars. 



IMFs, as for Case C. Accordingly, the IMFs that can 
reproduce the both CEMP star statistics fall in the nar- 
rowest range of M mc i = 2.5 — 7.0 Mq with rather small 
upper mass limit, almost irrespective of the dispersion, 
on the M md - A M diagram in FigO The CEMP-s star 
faction remains smaller than ~ 20% because of the con- 
tribution of the stars born as single. 

In conclusion, the statistics of CEMP stars demand 
the IMFs for the EMP population, peaking at the 
intermediate-mass stars or the massive stars, by far 
higher mass than those of Pop. I and II stars, irrespec- 
tively of the assumed mass-ratio distribution. The pres- 
ence of CEMP-nos stars in a significant number of the 
CEMP-s stars excludes the IMFs of small mass. The 
derived mass range varies by a factor of ~ 2, from the 
highest M m d > 7.1 Mq for the mass-ratio distribution of 
increase function of the mass ratio (Case B) to the lowest 
2.5 < M m( j/ Mq < 7 for the mass-ratio distribution of 
"independent" coupling (Case D). This tendency is ex- 
plained in terms of the difference in the averaged mass 
of the primary companion of the EMP survivors; if the 
contributions to the EMP survivors decrease rapidly with 
the mass of primaries, a relatively higher-mass IMFs re- 
sult without an upper mass limit imposed, while if the 
contribution to the EMP survivors are weakly dependent 
or independent on the primary masses, an upper limit is 
set with the relatively smaller-mass on the IMFs. 

2.3. Constraints from Galactic Chemical Evolution 

In this section, we discuss that additional constraints 
can be derived from the relationship between the total 
number of EMP survivors and the iron yields from the 



EMP population. 

Figure 0] shows the contour maps of the total stellar 
mass, Memp, of EMP population, in cq. (|10p . neces- 
sary to leave the observed number of EMP survivors on 
the M m d-AM diagram. The total stellar mass increases 
for higher mass IMFs to produce the given number of 
low-mass survivors. For the flat mass-ratio distribution 
(Case A; left panel), the total stellar mass of EMP pop- 
ulation is essentially determined by M m d in proportion 
to M^ d and only weakly dependent on Am for large 
M m d, since almost all EMP survivors formed as sec- 
ondary (i.e. second term of eq. [11] is dominant) for IMF 
of log(M md /0.8M Q ) > A M - Figure [5] shows the con- 
tours of the total iron mass, Mpe.EMP, produced by the 
massive stars of EMP population, in eq. (fT3|) . The iron 
production also increases with larger M m d but is more 
sensitive to Am since the dependence of the supcrnovae 
fraction differs across the border of M m d — M up ; the 
amount of produced iron increases (or decreases) with 
A M for given A/ md < M up (or M md > M up ), little de- 
pendent on Am for M m d — M up . 

For the "independent" coupling (Case D; left panel), 
the fractions of giant EMP survivor and supernova 
among EMP population stars are given by, 

/ G = (1 + / 6 )£(0.8)AM G (20) 
/sn = (l + fb) I dmi£(roi) (21) 

JM up 

The total mass and the amount of iron production of 
EMP population are sensitive both to M m d and Am, 
especially for small Am and large M m d in contrast to 
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Fig. 4. — The contour of the total mass of EMP population stars, Memp = 10 7 , 10 s and 10 9 Mq on the diagram of parameters (medium 
mass M md and dispersion Am) in the log- normal form for the fiat distribution (Case A: left pannel) of the mass ratio and independent 
counling (Case D: right pannel). Shaded area denotes the parameter range consistent with the statistics of CEMP stars. 

with the other cases. For small Am, therefore, the frac- 
tion of low-mass stars varies greatly with Mmd, and the 
both contours of Memp and Mf^Memp) converge to 
M m d — 2 — 3M Q . As Am increases, the differences 
from Case A diminish since the IMFs tend to extend 
into the low-mass stars, and in particular, for Am ^ 0.4 
and M m a < 3M Ql the contours in the both panels re- 
semble each other to run through the similar parameter 
spaces. 

From the comparison with the total amount of 
iron A/fc,emp, necessary for the chemical evolu- 
tion, in this diagram, the parameter space where 



M Fe ,EMp(M md , A M ) > M Fe ,haio = 10 5 5 M Q is excluded 
by the overproduction of iron or by the underproduc- 
tion of EMP survivors. For the parameter space where 
M Fei EMp(M md , A M ) < M FCiha i , on the other hand, the 
stars of EMP population can leave the number of EMP 
survivors currently observed but are short of iron pro- 
duction, so that the chemical evolution demands other 
sources of iron production without producing the low- 
mass stars that survive to date. For a flat mass-ratio dis- 
tribution, the IMFs that can satisfy the condition of iron 
production coincide the IMFs, derived above from the 
statistics of CEMP stars (shaded area) in the parameter 
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Fig. 5. — Constraints on the IMF of EMP population, derived from the number of EMP survivors and the iron production by the stars 
of EMP population on the Af md — Ajyj plane. Solid lines denote, the loci of IMFs which can produce the amounts of metal production, 
Mp„ emp = 10 4 ; 10 5 , and 10 6 Mq. Dashed lines denote the contour of carbon production by AGB stars to the iron production, [C/Fe] = 2, 
1, and 0. 



range of M md ~ 10 - 16 M© and A M ^ 0.3 - 0.6. For the 
"independent" coupling, the parameter range of IMFs 
that satisfy the condition of iron production also over- 
lap the shaded area of parameter range, derived above 
from the statistics of CEMP stars, but with the mass 
M mc i — 3.5 — 5.1 Mq, slightly smaller than for Case A and 
only for a small dispersion of Am — 0.23—0.35. For larger 
Am, even the highest-mass IMFs of M m d = 7.5 M© re- 
sult to be slightly short of, or marginally sufficient at the 
most, iron production. 

For the two other mass-ratio distributions of n oc q 
(Case B) and n oc 1/q (Case C), the iron production, 
Mpe,EMP, with a given IMF results to be larger or smaller 
than for Case A because of the difference in the number 
of massive stars exploded as supernova per low-mass sur- 
vivor (e.g., by a factors of 1.38 and 0.47, respectively, per 
a star of m = 0.8 M© and the IMF of M md = 10 Mq 
and Am = 0.4). The iron production then demands 
smaller-mass (higher-mass) IMFs for Case B (Case C) 
as compared with Case A, the shift of IMFs in an op- 
posite direction, discussed from the statistics of CEMP 
stars. Accordingly, for these two extreme cases, the pa- 
rameter ranges for the IMFs derived from the statistics 
of CEMP stars and the iron production are marginally 
overlapped (Case C) or are dislocated with a narrow gap 
(Case B), although a definite conclusion waits for future 
observations, in particular, to improve the estimate of 
total numbers of EMP stars (see Appendix). 

The relative production rate of carbon to iron may also 
impose additional constraint since the intermediate-mass 
stars enrich intcrgalactic matter with carbon through th e 
mass loss on the AGB, as discussed bv lAbia et al.1 (|2001h . 
In particular, when Am is small and M m( j is in the range 
of intermediate- and low-masses, the intermediate-mass 
stars much surpass the massive stars in number and eject 
more carbon than the latter eject iron. We compute the 
amount of carbon ejected by AGB stars by taking the 
carbon abundance in the wind ejecta of AGB stars at 
[C/H] = 0, and the remnant mass at 1M©. Contours of 
[C/Fe] = 2, 1,0 are plotted in the figure (dashed lines), 
for which only the carbon from the AGB stars are taken 
into account. The overabundance of carbon excludes the 
IMFs with low dispersion and low medium mass; it ex- 
cludes the parameter space in the range of Am < 0.2, 
derived by the CEMP star statistics for Case D, but 



has nothing to do with the high-mass IMFs derived for 
Case A. 

We demonstrate that the IMFs, derived from the ob- 
served properties of CEMP stars, have the parameter 
ranges that can explain the chemical evolution and the 
production of low-mass stars, consistent with the obser- 
vations, both for the flat mass-ratio distribution and for 
the "independent" coupling. In Appendix we will discuss 
the converse to demonstrate that the argument based on 
the total number of EMP survivors and the total iron 
production can potentially provides more stringent con- 
straint on IMFs, independent of the argument based on 
the CEMP star statistics. 

Relative abundances of other elements may also be 
affected by the IMF. Theoretical study of super- 
nova nucleosynthesis suggests the peculiarities and vari- 
ations of yi elds for the metal-free and extremely metal- 
poor s tars (IWooselv fc Weaver 119951; lUmeda fc Nomotol 
20021: iHeeer fe Wooslevi 12002 iTominaga et al.1 120071 : 



Heeer fc Wooslevil2008ft . The supernova yields are, how- 



ever, sensitive to the assumption of model parameters 
such as the explosion energy and the treatment of non- 
axisymmetric effects, and currently subject to the large 
uncertainties. In this paper, therefore, we are con- 
cerned with the iron yields as an indicator of the chemi- 
cal evolution, and defer detail study about the abundance 
pattern of various elements in future works. As for the 
iron yield, recently, type la supernovae with short de- 
lay time and their contributio n of iron production are 
discu ssed by some authors (e.g. IScannapieco fc Bildstenl 
|2005[) . although the evolutionary scenario is not yet clear. 
The iron reduction is suggested to be several times larger 
than by type II supernovae, and yet, will hardly affect 
our results since the ratio, Jsn/ fg depends strongly on 
the IMF. 



2.4. Distinctive Features of EMP Survivors 

The different assumptions on the mass-ratio distribu- 
tions admit the parameter ranges of high-mass IMFs that 
can reproduce the statistics of CEMP stars and the chem- 
ical evolution, consistent with the existent observations. 
The predicted mass ranges differ by a factor of 2 or more 
between M m( j — 5 — 20 Mq. Although hardly distinguish- 
able from the observations discussed so far, they surely 
make the differences in the properties of EMP survivors. 
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Fig. 6. — The mass function of EMP survivors under the different 
assumptions on the coupling mass distribution of binaries. Here the 
parameters of IMFs are taken to be M m ^ = 10 Mq and Am = 0.4 
(Cases A-C). 

We discuss the imprints that the mass-ratio distributions 
have left on the current EMP survivors and investigate 
the possibility of discriminating the mass coupling of bi- 
nary systems in the EMP population, especially for the 
two distinct distributions of the flat mass-ratio distribu- 
tion and the "independent" coupling. 

Firstly, an obvious difference is the mass distribution 
function of EMP survivors. For a given IMF, £(m), 
the mass distribution, £emp-sui-v (to) , of EMP survivors 
is given by; 

r-ra 

&MP-surv(»n) = (1 - fb)€(™>) + M(m) / n(m 2 /m)/mdm 



+ /b / £(mi)n(m/mi)/midmi. 

/O.8M 



(22) 



Here a low-mass binary, whose components are both less 
massive than 0.8 Mq, is counted as one object with the 
primary star. Figure [5] shows the mass distributions of 
EMP survivors (to < 0.8 Mq) for different assumptions 
of mass-ratio distributions Cases A-C. For these mass- 
ratio functions, the mass distribution of EMP survivors 
is nearly proportional to the mass-ratio distribution n(q) 
because almost all of them come from the secondary 
stars; the contribution from the primary components are 
denoted by thin solid line, and the same contribution 
comes from the stars born as single. For the "indepen- 
dent" coupling, in contrast, the ^EMP-surv(m), has the 
same form as the IMF and the number of EMP survivors 
decreases rapidly as the stellar mass decreases. 

Secondly, the fraction of double-lined binary and the 
contribution of stars born as single among EMP sur- 
vivors may differ according to the mass-ratio distribu- 
tion. The EMP survivors born as binaries are divided 
into three categories according to the mass of the pri- 
mary stars: (1) the low-mass binaries with the primary 
of mass mi < 0.8 Mq, (2) the white-dwarf binaries of 
primary stars of mass between 0.8 Mq < toi < M up , 
and (3) the supernova binaries of primary stars of mass 
mi > M up . The fraction of low-mass binaries with the 
primary stars of mass m < 0.8 Mq in the EMP survivors 
of mass between to to to + dm is given by 



XMB 



flat mass-ratio distribution, this gives a significant frac- 
tion of <y9 S m V ,LMB(0.8 Mq) = 7.3% for M md = 10 Mq and 
Am = 0.4, and increases with Am to 16% for Am = 0.5 
and with decreasing M m d to 18% for M mc j = 5M Q , re- 
spectively. We note that these values depend weakly on 
fb since most of the EMP survivors are from the bina- 
ries. The number of low-mass binary decreases rapidly 
for smaller masses while the number of EMP survivors, 
formed as the low-mass members of white dwarf binaries 
or supernova binaries, remains constant. 

For the "independent" coupling, the fraction of low- 
mass binaries in the EMP survivors reduces to; 



(m) = fb{£(rn)/rn] / n(rri2 / rn)drn2 / 



,EMP- 



,(m) 



They can be detected as double-lined binary. 



(23) 
For the 



,lmb m 



2fb 



£(TO 2 )dm 2 



(1 + fb) 



-2/ 6 



0.8 A/q 



£(mi)dmi 



(24) 



which gives a 

^surv,LMB(0.8 -Mq) 



much smaller fraction of 
1.6% for M md = 5M q and 
Am = 0.4 as compared with the flat mass-ratio distri- 
bution. The fraction may increase for smaller medium 
mass, to 5.5% at M m & = 3M©, and for larger dispersion, 
to 3.9 % and 9.5% at Am = 0.5 and 0.7, respectively, 
although these may cause underproduction of iron, 
in particular for smaller M m( j, as seen from Fig. [3] 
(bottom panel). In this case, the proportion of the EMP 
survivors, born as single stars, is fairly large as given by 



9: 



surv.sm; 



,(m) ~ (i-/ 6 ) 



(l + / 6 )-2/ h 



0.8 Mq 



^{mi)dm\ 
(25) 



Consequently, nearly one third of EMP stars were born 
as single stars, for /{, = 0.5, which is much larger fraction 
than in the case of the flat mass-ratio distribution. 

Thirdly, the fraction, (^surv.sNB, of supernova binaries 
with the primary stars of mass toi > M up also differs 
between the two mass-ratio distributions. For the flat 
mass-ratio distribution, almost all of the EMP survivors 
belong, or have been belonged, to the binary systems, 
and the fraction is given by 



^surv,sNB(™) = fb / n(n/mi)t(mi)/midmi/£ 



EMP- 



r(m), 



(26) 

and amounts to ~ 50%. For the " independent" coupling, 
on the other hand, one third of EMP survivors are single 
stars from their birth, and the percentage of supcrnovae 
binaries is relatively small, as given by 



y S urv,SNB(w) = 2/ 6 / £(mi)dmi / + 

>M up 

r0.8M e 

-2 / fc / £(mi)dmi 



(27) 



and turns out to be ~ 20%. The EMP survivors from the 
supernova binaries have experienced a supernova explo- 
sion of the erstwhile primary stars at close distances and 
are thought to suffer from some abundance anomalies, 
affected by supernova ejecta. Accordingly, these stars, 
in particular from the binaries of sufficiently small sep- 
arations, may be discriminated by a large enhancement 
of elements, characteristic to the supernova yields. 
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These differences in the properties of remnant EMP 
survivors may potentially serve as tools to inquire into 
the nature of EMP binaries and to distinguish the 
mass-ratio distributions. Among the EMP stars, sev- 
eral double-lined spectroscopic binaries are reported 
in the literature. If we restricted to the mctallic- 
ity range of [Fe/H] < —3, for which the observa- 
tions with high-resolution spectroscopy may be re- 
garded as unbiased, there are two stars CS22876-032 
([Fe/H] ~ -3. 6, V = 12.84, P = 424.7 d, and 
7Ti9 / mi ~ 0.89; iThorburn fe Beers"! 119931 ; iNorris et all 
[2000HGonzalez et al.H2008ft and CS 22873-139 ([ Fe/H] ~ 
-3.4 , V = 13.8, P = 19.165 d, m 2 /m 1 ~ 0.92; iPrestonl 
11994 12000t ISpite et al.ll2000D with the detailed analyses 
and one star HE 13 53-2735 ([Fe/H] ~ -3.2, V = 14.7; 
iDepagne et al.|[2000T ) without the binary parameter. So 
far 39 dwarf stars of [Fc/H] < —3 are confirmed by 
the high-resolution spectroscopy (we define the dwarf as 
log g[ cm s -2 ] > 3.5), and hence, the fraction of low-mass 
binaries, composed of two unevolved EMP stars, turns 
out to be > 3/39 ~ 7.7%. It seems to be compatible with 
the estimated fraction </?sui-v,lmb(0.8 Mq) — 7.3% for the 
frat mass-ratio distribution but to be a little larger, or 
marginal, for the "independent" coupling. We note, how- 
ever, that only with two samples, the above fraction may 
be subject to significant observational selection effects. 
These stars have to be concentrated near to the upper- 
end of main-sequence since they are found among the 
candidates, selected from the flux limited surveys, and 
the mass ratio has to be sufficiently large for the lines 
of two components to be observed. The actual fraction 
has to be larger than observed if we take into account 
the detection probability due to the orbital phase and 
to the inclination angle, and the rather narrow range of 
mass-ratios for the observed double-lined binaries. On 
the contrary, the larger survey volume by a factor up 
to 2 3 / 2 for the double-lined binaries due to the sum of 
luminosities may reduce the actual fraction. More obser- 
vations for the main-sequence EMP binaries and the bias 
corrections are necessary to discriminate the mass ratio 
distributions. 

It may be more straightforward to compare our results 
with the mass distribution function of EMP survivors. 
From the existent observations, however, it is rather 
hard to determine since the observed dwarfs are mostly 
concentrated near to the upper end of main sequence. 
An e xception is a carbon dwarf G77-61 of [Fc/H] = 
4.03 dPlez fe Cohenl l2005h whose mass is inferred at 
0.3 — 0.5 Mq, b ut it was found am ong the proper- motion- 
parallax stars (jDahn et al.lll977t ). not from the surveys. 
We have to wait for the larger-scaled surveys in near fu- 
ture to reveal the distribution of EMP survivors of low 
masses. As for the supernovae binaries, they are ex- 
pected to be related to the large star-to-star variations 
in the surface elemental abundances, in particular, with 
those of r-process elements, ranging more than by two 
orders of magnitude. It is necessary, however, to un- 
derstand the nature of interactions between the super- 
nova ejecta colliding at very high velocity and the near- 
by low-mass stars before the meaningful conclusions can 
be drawn from the observations. 

3. METALLICITY DISTRIBUTION FUNCTION OF EMP 
STARS 



We have shown that the high-mass IMFs with the bi- 
nary provide a reasonable explanation of the observed 
properties of EMP stars in the Galactic halo, revealed 
by the recent large-scaled HK and HES surveys. In 
this section we discuss the consequence of derived IMF 
on the metal enrichment history of Galactic halo up to 
[Fe/H] = —2.5 to study their relevance to the metallic- 
ity distribution function (MDF), observed for the EMP 
stars. As seen above, the derived IMFs using the dif- 
ferent assumptions on the mass-ratio distribution func- 
tion are hardly distinguishable by the current observa- 
tions. In the following, therefore, we assume the IMF 
with the medium mass M m d = IOMq and the disper- 
sion Am = 0.4 in the lognormal form, which is at least 
compatible with any of the assumed mass-ratio distribu- 
tions. 

3.1. Simple Model of Chemical Evolution 

Under the assumption that matter ejected from su- 
pernovae spreads homogeneously and is recycled instan- 
taneously, the iron abundance, X Fe , of our Galaxy of 
(baryonic) mass Mh can simply be related to the cu- 
mulative number, N(Xp e ), of the stars born before the 
mctallicity reaches X Fe as; 



M h X Fc = {Y Fe )N(X Fe )f SN . 



(28) 



where (Y Fe ) is the averaged iron yield per supernova and 
/sn is the fraction of EMP stars that have exploded as su- 
pernovae, defined in eq. (|14[) . By differentiating it with 
respect to [Fe/H] = \og(X FB /X Fe: Q), the number dis- 
tribution of EMP survivors is written as a function of 
mctallicity in the form 



n([Fe/H] 



dN{X_ 



Fc I 



d[Fe/H] (F F c)/sn 



ln(10)A FcQ 10 [Fc/H1 . 

(29) 

This shows that the number distribution of EMP sur- 
vivors is simply proportional to the iron abundance apart 
from the variation of (Yp e ) through the IMF and the lat- 
ter is s mall enough to be neglected for M m <j < 20 Mq 
(see Fig. lAll in Appendix). 

Figure [7] depicts the number distribution of EMP sur- 
vivors and compares it wi th the observed M DF pro- 
vided by the HES survey (jBeers et alJ 12005! ). We as- 
sume stars of mass m > M up = 8 M Q become type 
II supernova and eject (Y Fc ) = 0.07 Mq of iron. In 
this figure, the theoretical MDF, v SUIV , is evaluated 
under the same flux-limited condition as the observed 
MDF is derived; i/ surv ([Fe/H]) = n([Fe/H])/ G x (40%) x 
(8225 degree 2 /An sr) x 1.93. Here the fraction of follow- 
up observation and the sky coverage are taken into ac- 
count: as for the contribution of TO stars, we take the 
same ratio to the giants as in the observed sample under 
the assumption that the giant survivors are all reached 
in the survey area. Solid line shows the MDF for the 
IMF with M md = 10M o and A M = 0.4 with the 50% 
binary fraction, derived above for EMP population stars 
for the flat mass-ratio distribution, and it is similar to 
the other mass-ratio functions, as discussed in §2.3. This 
reasonably reproduces the observed MDF between the 
metallicity —4 < [Fc/H] < —2.5, as expected from the 
discussion in the previous section. 

In this figure, we also plot the MDF using the low- 
mass IMFs, the Salpeter's power-law mass-function as 
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M md =10,A M =0.4 
Salpeter 

10000 - Lucatello{2005) 




Fig. 7. — Comparison of the theoretically predicted MDF from 
the IMF derived from the statistics of EMP stars (solid line) with 
the observed MD F obtained by the HES survey (shaded columns; 
Beers ct al. 2005). The number distribution, u Burv is computed 
with the log-normal IMF of medium mass M mt j = 10 Mq and dis- 
persion Am = 0.4 with the 50% binary fraction under the same 
flux-limited condition as the HES survey. The examples of low- 
mass IMFs are shown for the Salpeter's power-law mass function 
(dashed line) and for the log-normal IMF (M mr) = , 79 M p; and 
Am = 1.18, dotted line), derived by[Lucatcllo et al. (2005]) from 
the CEMP-s star statistics alone, which bring about the overpro- 
duction of low-mass survivors. 

observed among the present-day stellar populations and 
th at derived only f rom t he statistics of CEMP-s stars 
bv iLucatello et all (|2005l M md = 0.79 M© and A M = 
1.18/ In 10). They bring about the overproduction of 
EMP survivors by a factor of more than a few hundreds 
not only from the low-mass members of binaries but also 
from the primary stars and the single stars; both the 
IMFs give the similar MDF since our flux-limited sam- 
ples are dominated by the giants and luminous dwarfs of 
mass M ~ 0.8 Mq. This means that the EMP survivors 
is by far a small population as compared with the stellar 
systems of Pop. I and II, and it is only with the high- 
mass IMFs that can make the EMP population produce 
sufficient amount of metals to enrich the early Galactic 
halo without leaving too many low-mass survivors now 
observable in Galactic halo. 

In addition, we see in this figure that the slope of ob- 
served MDF is consistent with the prediction from the 
simple one-zone approximation at least for [Fe/H] > —4. 
It implies that the IMFs have little changed while the 
Galactic halo has evolved through these metallicities. 
Beyond [Fe/H] ~ —2, the observed MDF derived from 
the HK and HES surveys seems to be underestimated 
since those objects are out of the metallicity range sought 
after by the survey and subject to imperfect selection. 

3.2. Effect of Hierarchical Galaxy Formation 

The observed MDF of Galactic halo stars has a sudden 
drop at [Fe/H] < —4, and only three stars are found 
below it 3 . We propose the mechanism responsible for this 
depression of low- metallicity stars from the consideration 
of the Galaxy formation process. 



3 There are two more stars with the iron abundances reported 
below [Fe/H] < -4 ; CD -38°245 with [Fe/H] = -4.19 ± 0.10 
I Cavrel et atl |2004|) and G77-61 with [Fe/H] = -4.03 ± 0.1 
i PIez fc Cohe n 2005). We will omit these two stars in our 
discussion since larg er abundances of [ Fe/H] = —4.07 ± 0.15 
lIFrancois et al J 120031) and -3.98 ± 0.15 l INorris et al.|[200ll) have 
been reported for the former, and hence, their abundances are 
closer to the EMP stars of [Fe/H] > —4 than to the other three 
HMP/UMP stars. 



In the current cold dark matter (CMD) model, galax- 
ies were formed hierarchically They started from low 
mass structures and grew in mass through merging and 
accreting matter, finally to be large-scale structures like 
our Galaxy. In the hierarchical structure formation sce- 
nario with ACDM cosmology, the typical mass of first 
star forming halos is ~ 1O 6 M for the dark matter and 
~2x 10 5 Mq for t he gas (e.g., see lTegmark et alJll997t 
iSpergel et alTl2007h . 

In these first collapsed gas clouds, the first stars con- 
tain no pristine metals except for lithium. When the first 
star explodes as supernova, it ejects ~ 0.07 Mq of iron, 
which enriches the gas cloud of mass ~ 2 x 10 5 Mq where 
it was born up to the metallicity of [Fc/H] ~ —3.5 if the 
ejecta is well mixed in the gas cloud. We call this event 
the "first pollution". Consequently, the 2nd generation 
stars have the metallicity of [Fe/H] ~ —3.5. 

In the course of time, the mini-halos that host the gas 
clouds merge with each other and accrete the intergalac- 
tic gas to form early Galactic halo with the baryonic 
mass of 10 11 M Q . We may take the metallicity of this 
early Galactic halo to be [Fc/H] ~ —4 because of the 
scarcity of stars of metallicity [Fe/H] < —4. The cu- 
mulative number of stars born before the early Galac- 
tic halo is enriched up to [Fc/H] = —4 is estimated at 
JV(lO _4 AF e 0) = 7 x 10 5 with taking into account the su- 
pernova fraction /sn • If the mini-halos of larger masses 
stand between the first collapsed halos and the Galac- 
tic halo, the dilution of iron with unpolluted primordial 
gas can give birth to the stars of smaller metallicity of 
[Fe/H] ~ —4, and then, the metallicity at the formation 
of Galactic halo can be larger to increase the cumulative 
number of stars in accordance (see below). 

We may estimate the fractions of both the first gener- 
ation stars without metals and the 2nd generation stars 
of the metallicity [Fe/H] ~ —3.5, respectively, assuming 
that stars are born with an equal probability whether in 
the gas clouds, polluted with metals, or in the primordial 
gas clouds. Here it is worth noting that some recent 
computations of star formation demonstrate that the 
low-mass stars can be f ormed as the binary members even 
out o f metal-free gas (|Clark et al.l 12008: Machida et alj 
|2008| ). Accumulated number, Np op m, of Pop III stars, 
born of gas unpolluted by SN ejecta, when the average 
metallicity reaches Xp c , is given by 



Am 



M h 



PopIII 



M C / SN 



1 — cxp 



M r X 



Fc 



(*Fe> 



(30) 



where M c (= 2 x 10 5 Mq) is the mass of gas in the first 
star forming clouds. If we assume the same IMF and 
binary parameters as in the stars of EMP population, 
then, we expect that the number of Pop III stars is 

AWm^O^Fc,©) = 3.1 x 10 5 , (31) 
and the number of Pop III survivors is 

,-0.8 Mq r 

3.1 x 10 5 x / dm[£(m) + /& / n(m/mi)£(mi) 



0.08 Mq 



0.8 M e 



= 1.3 X 10 4 , 

and similarly we have 5.5 x 10 4 and 2.3 x 10 3 of the 2nd 
generation stars and their survivors, formed before the 
averaged metallicity of the Galaxy reaches [Fe/H] = —4. 
The IMF of Pop. Ill stars may differ from EMP stars 
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but the existence of the stars with [Fc/H] < —5 suggests 
that the low-mass stars can be formed before the first 
pollution. 

Figure [H] illustrates an expected MDF with the hierar- 
chical structure formation. After the formation of large 
Galactic halo, the metal enrichment process is thought 
to follow the argument of the previous subsection. Thus, 
we can explain the cutoff around [Fc/H] ~ —4 natu- 
rally. Shaded columns indicate the initial distributions of 
Pop. Ill stars and of the 2nd generation stars formed in 
the low-mass clouds. The 2nd stars were mixed and ob- 
servationally lost their identities among the stars formed 
in the merged halo. On the other hand, Pop. Ill stars 
should form a distinctive class. From the above estimate, 
we expect ~ 23 Pop. Ill survivors in the existing flux- 
limited samples of HES surveys. It is true, however, that 
there is no star with zero metallicity among the stars 
detected by the existent surveys. We may propose one 
scenario to explain this absence that Pop. Ill survivors 
are no longer remain metal-free at present since their sur- 
face are polluted by the accretion of interstellar matter, 
enriched with metals ejected by the supernovae of the 
first and subsequent generations. With the surface pol- 
lution of [Fc/H] ~ —5, they are observed as HMP stars. 
We discuss about evolution of Pop. Ill stars with pollu- 
tion in £)4.1I Similarly, the number of second-generation 
of stars is estimated at ~ 4 in the same flux-limited HES 
sample, indicative that most of EMP stars are formed 
of mixture of the ejecta from plural supernovae. This 
has direct relevance to the study of the nucleosynthetic 
signatures on the EMP survivors and the imprints of su- 
pernovae of the first and subsequent generations. 

4. CONCLUSIONS AND DISCUSSION 

We have studied the initial mass function (IMF) and 
the low-mass star formation with the chemical evolution 
of the Galactic halo population on the basis of the charac- 
teristics of extremely metal-poor (EMP) stars, revealed 
by the recent large-scaled HK and HES surveys; the ob- 
servational facts that we make use of are; (1) the over- 
abundance of carbon-enhanced EMP (CEMP) stars, (2) 
the relative frequencies of CEMP stars with and without 
the enrichment of s-proccss elements, (3) the estimate 
of surface density or total number of EMP stars in the 
Galactic halo, and (4) the metallicity distribution func- 
tion (MDF). We take into account the contribution of 
binary stars properly, as expected from the younger pop- 
ulations. In Paper I, the high mass IMF peaking around 
~ 10 Mq is derived for the stars of EMP population and 
it is shown that the binary population plays a major role 
in producing the low-mass stars that survive to date, 
but by using the flat mass-ratio distribution between 
the component stars. In this paper, we examine these 
properties of the stars of EMP population and EMP sur- 
vivors for the different types of mass-ratio distributions 
and investigate the constraints on the IMFs of the stars 
of EMP population and discuss the observational tests 
of discriminating them. The derived IMFs are applied to 
understand the characteristics of MDF and the nature of 
EMP stars including HMP/UMP stars, provided by the 
surveys. 

Our main conclusions arc summarized as follows; 
(1) The statistics of CEMP stars are explained by the 
high-mass IMFs with the binaries of significant fraction. 



Predicted typical mass is significantly larger than Pop- 
ulation I or II stars, irrespective of the assumptions of 
the mass-ratio distribution. The mass-ratio distribu- 
tion with a preference for nearby equal masses demands 
the IMF with higher typical mass M m( j > 7 Mq and 
smaller dispersion (Am < 0.6). While the mass-ratio 
distribution in favor of smaller mass secondary or the in- 
dependent combination of two stars with the same IMF 
demands smaller typical masses around M m d ~ 5M 
irrespective of Am- 

(2) High mass IMFs with M md ~ 5 - 20 Mq derived 
from the statistics of CEMP stars agree with those de- 
rived from the low-mass star formation and the chemical 
evolution of Galactic halo based on the number of gi- 
ant EMP sourvivors evaluated from the surveys. IMFs 
with M m( j 3> 20 Mq are excluded by overproduction of 
iron or underproduction of EMP survivors. IMFs with 
M m d -c 5 Mq need other iron source(s) without produc- 
ing low-mass stars. 

(3) The mass-ratio distribution of binaries in the EMP 
population can be discriminated by the imprints left on 
the EMP survivors such as the mass function, the binary 
fraction, and the fraction of stars influenced by the su- 
pernova explosion of primary stars. In particular, the 
flat mass-ratio distributions predict significant fraction 
(~ 7.7%forM m a = 10) of double-lined spectroscopic bi- 
naries while the mass-ratio distribution of "independent" 
coupling predict much lower fraction. Among the 39 un- 
evolved stars of [Fe/H] < —3, studied spectroscopically 
to date, three double-lined binaries are found, but there 
may be significant uncertainties and biases for the ex- 
istent surveys and future observations can discriminate 
the distributions. 

(4) The observed MDF of EMP survivors is consequent 
upon the derived IMF with the contribution of the bina- 
ries. There is no indication of significant change in the 
IMFs between the metallicity of -4 < [Fe/H] < -2. The 
depression of stars below [Fe/H] < —4 is naturally expli- 
cable within the current favored framework of the hier- 
archical structure formation model. Then, the Pop. Ill 
stars born of primordial gas, and also, the stars in the 
primordial clouds before they are contaminated by their 
own supernovae, should form the distinct class other than 
EMP stars, and may have the relevance to HMP and 
UMP stars observed at lower metallicity, as discussed 
below in this subsection. 

The feature of our approach is to take into account 
the stars born in binary systems properly in discussing 
the low-mass star formation in early Universe, based on 
the finding in Paper I. In addition, we make full use 
of available information from the existent large-scaled 
surveys and to draw the maximal constraint on the 
early evolution of our Galactic halo. The known EMP 
stars ([Fe/H] < —2.5) with the detailed stellar param- 
eters amount to ~ 400 in number (SAGA Database; 
ISuda et aTll2008[ ). and allow us to discuss the averaged 
properties as studied in this paper. The existent sur- 
veys reach sufficiently deep and the nominal depth in 
the magnitude 12 < B < 17.5 may corresponds to the 
heliocentric distance d ~ 10 kpc or beyond for giants of 
L ~ IOOLq while d ~ 0.3 — 3 kpc for dwarfs. Incomplete 
coverage may affect the estimate of the total number of 
EMP survivors in the Galactic halo, and yet, the uncer- 
tainties will not be so large to exceed an order of mag- 
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nitudes, judging fro m the density d istribution of bright 
halo p ~ r~ 3r 3 5 (|Majewskilll993f l. Furthermore, be- 
cause of strong dependence of the number ratio between 
the low-mass stars that survive to date to the massive 
stars that have exploded as supernovae on IMF, our dis- 
cussion in §2.3 through the iron production consistent 
with the number of EMP survivors will be left largely 
unaffected even quantitatively. In order to improve and 
sharpen our conclusions, we have to wait fo r the future 
large r-scaled surveys s uch as SPSS /SEG UE (jBeers et alJ 
l200l and LAMOST (|Zao et al.llBM) . Also the high 
dispersion spectroscopy is necessary to understand the 
characteristics of EMP stars. 

The constraints on the IMFs derived in this work may 
serve as the basis of understanding the formation and 
early evolution of the Galaxy. Although the HK and 
HES surveys may suffer from i mper fect selections for 
[Fe/H] > -2, iRvan fc Norrisl (|1991| ) report that the 
mctallicity distribution extends continuously from the 
peak at [Fe/H] ~ -1.6 down to [Fe/H] ~ -3 with- 
out a break for the halo subdwarfs in the kincmatically- 
selected samples. This may be taken in turn as suggest- 
ing that the high-mass IMF derived for the EMP popu- 
lation persists without indication of significant changes 
up to such a large mctallicity. It forms a contrast with 
the prevalence of low-mass IMF with the characteris- 
tic mass M ~ O.3M in the diverse conditions includ- 
ing the Galactic sphero id population o f the mctallicity 
[Fe/H] ~ -1.7 - -1.4 (|Chabriedl2003l) and the globu- 
lar cluster, while evidences, suggesting IMF biased to- 
ward the high masses , are reported fr om cosmological 
observations (e.g., see lElmegreenl [20081 ) . Accordingly, 
there should be the transition from the high-mass IMF 
to the low-mass one. Our result suggests that the tran- 
sition is postponed until high metallicity even beyond 
[Fc/H] ~ — 2 is reached. It is likely that the transi- 
tion may not be simply determined by the metallicity 
alone, and may be related to the structure changes as- 
sociated with the merging processes. For the proper un- 
derstandings of the transition, more theoretical works 
are necessary with the effects of hierarchical formation 
of Galactic structures properly taken into account. In 
discussing the primordial stars or HMP/UMP stars in 
the present work, we assume the metallicity at the for- 
mation of Galactic halo at [Fe/H] ~ —4. The detailed 
chemical evolution with the merger history taken into 
account is discussed in a subsequent paper (Komiya et 
al. 2008, in prepar ation), i n the similar w a ys as done 
bv lTumlinso"n1 (|2006l ) and bv lSalvadori et all |2007D . but 
taking into account the high-mass IMF, derived above, 
and the contribution of binaries. 

4.1. Origin of HMP/UMP Stars 

We end by discussing the consequences of the present 
study on the understanding of the origin of stars found 
below the cut-off of MDF. 

In our model, the stars made after the first pollution 
have the metallicity [Fe/H] ~ —3.5 and the stars with 

slightly lower metallicity of [Fe/H] ~ —3.5 4 are made 

in the merged clouds where metals are diluted with the 
primordial gas unpolluted by supernova ejecta. After the 
halos merge, the 2nd generation stars mingle and obser- 
vationally lose their identities among the stars formed in 
the merged halo. This means that the HMP and UMP 



stars are the stars formed before the first metal pollu- 
tion by the type II supernova in their host mini-halos. 
One possible scenario for these stars is that they are the 
survivors of Pop. Ill stars and the metal abundances at 
their surface arc influenced not only by the matter from 
their binary companions but also by the interstellar mat- 
ter, enriched with metals ejected by the supernovae after 
their birth. In fact, it is shown that their peculiar 
abundance patterns of light elements from Li, carbon 
through aluminum, including s-process elements such as 
Sc and Sr, observed for three known HMP/UMP stars 
can be reproduced by neutron-capture nucleosynthesis 
during the AGB phase of their p rimary stars even unde r 
the pristine metal-free condition (iNishimura et al.ll2008l) . 
As for iron group elements (Ca - Zn). ISuda et all (|2004h 
argue the effects of surface pollution of Pop. Ill stars 
through the accretion of interstellar gas to show that 
the main-sequence Pop. Ill stars can be polluted to be 
[Fc/H] ~ —3 while the giants to be [Fe/H] ~ —5 since 
the pollutant is diluted by the surface convection deep- 
ening ~ 100 times in mass on the giant branch. Thus, 
the Pop. Ill survivors have evolved to giants to be ob- 
served as HMP/UMP stars. As for a sub-dwarf HMP 
star HE1327-2623, the dilution of the accreted iron group 
elements has occurred in the envelope of primary star on 
the AGB because of the low-mass nature of its primary 
star (M < 1.5 M g . INishimura et aT1l2008h . 

Majority of the Pop. Ill survivors have also to be the 
secondary members of binary systems similar to the EMP 
survivors if their IMF and binary parameters are similar 
to EMP stars. Then some of Pop. Ill stars become 
carbon-enriched HMP/UMP stars with [Fe/H] ~ -5 
through binary mass transfer. If the mass of primary star 
is 0.8 M Q < mi < 3.5 M Q and 3.5 M Q <mi < M up , the 
primary star enhances the surface abundances of carbon 
and nitrogen though the He-FDDM and of carbon and/or 
nitrogen through TDU and hot bottom burning in the en- 
velope, respectively, which are transferred onto the sec- 
ondary stars through the wind accretion. It is to be noted 
that the primary stars of mi > 2 Mq have the accreted 
pollutants mixed inward into the whole hydrogen-rich 
envelope at the second dredge-up, and thereafter, evolve 
like the stars with the pristine metals. At the same time, 
the accreted matter is diluted in the envelope and the 
iron abundance is reduced to [Fc/H] ~ — 5 in the pri- 
mary stars. We estimate that ~ 35% of Pop. Ill stars 
become carbon-rich HMP /UMP stars under the same as- 
sumptions on the binary parameters as in Paper I. The 
surface abundances of main sequence stars can be smaller 
than stated in Paper I, however, since the accreted mat- 
ter mixes down and are reduced by an order of mag- 
nitude if the diffusion and thermohalin e mixing works 
(|Weiss et al.ll2000t [Stancliffe et al.ll200l . 

In Fig. [8j solid lines denote the expected MDF at 
the present days with the surface pollution taken into 
account. The basic form of observed MDF is repro- 
duced, i.e., the cutoff around [Fe/H] ~ —4, the scarcity 
of stars for the metallicity below it and the existence 
of a few HMP/UMP stars. From the above estimates, 
there should be ~ 23 Pop. Ill stars in the existent flux- 
limited samples of HES surveys; about a half of them 
may be discovered as giants with the surface metal pol- 
lution and one third as carbon stars. In actuality, only 
three HMP/UMP stars (two giants and one sub-dwarf) 
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Fig. 8. — Schematic drawing of MDF of EMP and Pop. Ill sur- 
vivors, constructed based on the hierarchical scenario for structure 
formation. Primordial main sequence stars with Z = polluted 
through the after-birth accretion of the interstellar gas, enriched 
with iron ejected by supernovae of the first and subsequent genera- 
tions, upto [Fe/H] ~ —3, and pollutants is diluted to [Fc/H] ~ —5 
in the surface convection of EMP survivors as it develops during 
ascent of red giant branch. 

are found to date, all enriched with carbon, among 153 
stars of [Fe/H] < —3, registered in SAGA database 
(iSuda et al.ll2008h . Since such low metal abundances can 
be discriminated only with high dispersion spectroscopy, 
~ 4.6 HMP/UMP stars are expected among the whole 
HES samples of 234 stars of [Fc/H] < -3. The observed 
numbers are significantly smaller than prediction from 
our model. The above estimates are made, however, un- 
der the assumption that the Pop. Ill stars are formed in 



the same IMF as EMP stars and with the same binary 
parameters. This may not be warranted and rather we 
may take that this deficiency may suggest a still higher- 
mass IMF and/or less efficiency of binary formation for 
Pop. Ill stars than the EMP stars. 

In the above discussion, we assume the closed box 
chemistry in the collapsed object before merging. It 
is shown that the hypcrnovae, exploded with a large 
energy of 10 52 erg, blow off the first co l lapsed objects 
of mass M ~ 10 6 M Q (jMachida et all 120051) ; if the 
first stars are sufficiently massive, the metal yields are 
spread into larger masses, and pollute the ambient gas 
be fore they collap s e to f orm mini-haloes, as discussed 
by ISalvadori et all (|2007f ). After that, the first stars in 
the collapsed clouds are no longer metal-free. Never- 
theless, those stars which are formed before each col- 
lapsed clouds are polluted by their own supernova form 
a distinct class from those which suffer from the first 
pollution. Further study is necessary to make clear the 
present appearance of the possible Pop III survivors and 
to settle the origin of HMP /UMP stars, in particular, for 
tiny amounts of iron-group metals and the overwhelming 
carbon-enhancement, shared by all these stars known to 
date. 

We benefit greatly from discussion with Dr. W. Aoki. 
This paper is supported in part by Grant-in- Aid for Sci- 
entific Research from Japan Society for the Promotion of 
Science (grant 18104003 and 18072001). 



APPENDIX 

LOW-MASS STAR FORMATION AND IRON PRODUCTION BY EMP POPULATION 

One of the important findings of the recent large-scaled surveys is the scarcity of EMP stars in the Galactic Halo. The 
HES survey gives the total number of EMP stars in our Galactic Halo at oemp — 796 sr" 1 (giants of ctemp.g — 412 sr" 1 
in eq. (jTTJ) plus turn-off stars <temp,to — 384 sr" 1 ) within the limiting magnitude B ~ 17.5. Similarly, the HK survey 
gives (Temp — 528 sr _1 within the limiting magnitude of B ~ 15.5; 114 stars of [Fc/H] < —3 are found by the medium- 
resolution, follow-up spectroscopy of 50% of the candidates, selected from the objective prism survey covering the 
2800 deg 2 and 4100 deg 2 areas in the Northern and Southern Hemisphere (|Beers fc Chris tlicb 2005). Because of the 
significantly large areas covered by these surveys (~ 20% of all sky with the follow-up observations), we may place 
reliance on these results, granted that they may not be complete. This also constrains on the IMF of stellar population 
that promoted the chemical evolution, or more specifically, the formation of metals and the low-mass survivors. In 
the paper, we have discussed the chemical evolution starting with the statistics of CEMP stars. In this Appendix, we 
show that the chemical evolution with the total number of EMP survivors provides more stringent constraints on the 
IMF of EMP population with the aid of the amount of ejecta from supernova models, independently of the statistics 
of CEMP stars. 

Our basic premise is that the same stellar population is responsible both for the production of metals and of low- 
mass survivors. In discussing the low-mass survivors, it is indispensable to take into account the contribution from the 
binaries. This is one of the major conclusions in Paper I. We assume that the stars are born not only as single stars 
but also as the members of binaries in an equal number and with the primary stars in the same IMF as the single 
stars. For a given IMF, then, the total number, AEMP.surv of EMP survivors, currently observed in the Galactic halo, 
is related to the cumulative number, -/Vemp, of stars of EMP population as; 

A^EMP.surv = A^EMp/surv 

/•o.8M pi dm-] 

= N EM p dm[£(m) + f b n{mlra x ) ], (Al) 

J Jo.SM e m i 

and hence, to the cumulative number of EMP supernovae as A^emp.sn = A'emp/sn = A r EMP,surv(/sN//surv)- These 
supernovae have to supply the amount of iron, Mf c ,emp in eq. (|18[) . in order to enrich the gas in the Galaxy of mass Mh 
with iron to promote the chemical evolution up to the metallicity [Fc/H] = —2.5. Then, we may derive the averaged 
iron yield, (IrJemp, per supernova of EMP population, necessary to explain the chemical evolution of Galaxy, by the 
relation (Yes)emp = - ^Fc.e mp /A^emp.sn for an assumed IMF with the mass-ratio distribution function. 

We show in Figure [All the averaged yield, (Yp e )EMP 5 as a function of M m d for Am = 0.4: upper panel for the 
observations of EMP stars of different evolutionary stages from the HES survey and of the total EMP stars from 
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Fig. Al. — The average iron yields per supernova, (lp e }EMPj demanded from the chemical evolution of Galactic halo that leaves the 
EMP survivors consistent with the observe d flux-limited samples, an d the th eoretical iron yields , (VjsOsNi c omputed from the theore tical 
supernova models with use of ion yields by Umcda & Nomoto (2002), and by Wooscly k Weaver (1995) and Hcgcr & Woosley (2002) as a 
function of M mt j with Am = 0.4. Top panel compares (yp e )EMP computed for the different samples of EMP stars, i.e., giants, dwarfs and 
total stars from the HES survey and total stars from HK survey, with the fiat mass ratio function, while the bottom panel compares the 
results with the different mass-ratio distributions for the giant samples of the HES survey. 

the HK survey with use of the IMFs with the flat mass-ratio function, and lower panel for the different mass-ratio 
functions with use of the observation of EMP giants from the HES survey. In order to compare the stars of different 
evolutionary stages, we include the effects of the limiting magnitude of the surveys by assuming the de Vaucouleurs 
density distribution, p oc exp(— r 1 / 4 ), with the radial distance, r, from the Galactic center, the same as the stars in the 
Galactic halo and by assigning the luminosity of L = Lq(M/ Mq) 3 - 5 and lOOi© to dwarfs and giants, respectively. 

The amount of iron demanded by the chemical evolution turns out to be a steep decrease function of M m( j since in 
order to leave a fixed number of low-mass survivors, the total number of stars of EMP populations, and hence, the 
supernova fraction increase rapidly with M m( j in particular near M m d ~ M up . The necessary yields computed from 
the different samples in upper panel show a fairly good agreement with each other. The difference between the giants 
and dwarfs for the HES samples is indicative of a relatively deficiency of dwarf stars compared with giants by a factor 
of ~ 2.2 in number, which may be attributed to rather crude assignment of averaged giant luminosity, and/or to the 
different efficiency of identifying giants and turn-off stars in the survey plates, and/or to the uncertainties in the spatial 
density distribution. The results for the HK survey and the HES survey also agree within the difference by a factor 
of ~ 2.3 in number despite the difference in the limiting flux by 2 mag, and hence, to the difference in the searched 
volume by a factor of ~ 20. 

The variations with the mass-ratio functions in the lower panel are caused by the difference in the number of super- 
novae per EMP survivor. As compared to the flat mass-ratio function, the mass-ratio function increasing (decreasing) 
with q give a larger (smaller) number of supcrnovae to produce one EMP survivors; the difference of which increases 
for higher-mass IMFs. 

These iron yields necessary to promote the chemical evolution may be compared with the theoretical iron yields 
predicted from the supernova models. The IMF-weighted iron yields, (iFe)sN, per supernova is given by using the iron 
mass, Ypcfju), ejected from a massive star of initial mass m as; 



SN 



Im up dmi£(mi)[Y Fe (mi) + f b J Mup/mi Y Fc (m 2 )n(q)dq] 
Im up dmi£(mi)[l + f b fif ar/mi n{q)dq] 



(A2) 



The IMF averaged yield (IrJsn is also shown in this figure, for which the theore tical yields are tak e n from the 
metal-deficient supernov a models computed by lUmeda fc Nomotol (|2002t ). and by IWooselv fc Weaver! (|1995f ) and 
iHeger fc Wooslevl ([2002]) . It is a slowly increase function of M m( j for M m d < 20 Mq with the increase in the fraction 
of more massive stars that ended as supernovae, while beyond it, the gradient grows steeper owing to the contribution 
of the electron pair-instability supernovae of M > 100 Mq. 

The averaged yields, demanded by the chemical evolution, and the theoretical IMF-weighted iron yields both meet 
with each other near M m( j ~ 6 — 11 Mq and with the iron yield (iFe) — 0.04 — 0.06 Mq per supernova. As typically 
seen for the flat mass-ratio distribution, the parameter range coincides with that we have derived for the IMFs from 
the CEMP statistics in Figs. [l][3] For higher-mass IMFs, the EMP stellar population cannot produce the sufficient 
number of low-mass survivors by themselves, while for lower-mass IMFs, it results short of iron production. The 
differences arising from the mass-ratio distributions seem discernible but not large enough to differentiate these mass- 
ratio distributions in view of the uncertainties of current observations. As compared with the flat distribution, the 
mass-ratio distribution increasing (decreasing) with q demands smaller (larger) M m( j, the opposite tendency derived 
from the CEMP statistics. These distributions prefer smaller (larger) number of EMP survivors, larger (smaller) 
fraction of CEMP-s stars and smaller (larger) ratio of CEMP-nos to CEMP-s stars. In principle, however, we can 
discriminate the mass-ratio functions in the EMP binaries, including those destructed already by the evolution, with 
use of the survey and observations of EMP stars in sufficiently large number and with sufficient accuracy, which waits 
for future works. 



18 



Komiya et al. 



In summary, the observed surface density of EMP stars indicates the high-mass IMF for the stars of early stage 
of Galactic evolution independently of the CEMP star statistics, and also indifferently of the assumed mass ratio 
distribution function. It is true that the current estimate of Nemp /, surv ma Y be subject to significant errors, and 
yet, this result is robust because of a strong dependence of (Ype on M m d, as seen from the figure, which gives 
<51ogA'/ m( j ~ 0.4 log TVemp surv In the above discussion, we assume a single log-normal IMF with the binary fraction 
for the stellar population. It is possible to assume the bi-modal IMF and to explain the production of iron and the 
formation of low-mass stars, separately, in terms of the combination of two stellar populations, one with a higher-mass 
IMF responsible for the iron production and the other with a lower-mass IMF for the low-mass survivors, respectively. 
In the case of bi-modal IMFs, the constraints, derived here, place an upper mass limit to the IMF of lower-mass 
population and an lower mass limit to the IMF of higher-mass population. It is to be noted that the IMF with the 
binary mass function of Cases A-C is regarded as a sort of bi-modal IMF with the primary plus single stars as the 
higher-mass population and the secondary stars as the lower-mass population (see Fig. 12 in Paper I); the separation 
of two IMFs differs with the mass-ratio function and the relative contributions of two populations vary with the binary 
fraction. In any case, as for the EMP stars in the Galactic halo, the statistics of CEMP stars, in particular, the ratio 
of CEMP-nos and CEMP-s stars, place the lower mass limit to the IMF of the lower-mass population, and hence, 
endorses the high-mass IMF, which narrows, if any, the contribution of the higher-mass population. 
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